Multiple electron trapping in the fragmentation of strongly driven molecules 
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We present a theoretical quasiclassical study of the formation, during Coulomb explosion, of two 
highly excited neutral H atoms (double H*) of strongly driven H2. In this process, after the laser 
field is turned off each electron occupies a Rydberg state of an H atom. We show that two-electron 
effects are important in order to correctly account for double H* formation. We find that the route 
to forming two H* atoms is similar to pathway B that was identified in Phys. Rev. A 85 011402 
(R) as one of the two routes leading to single H* formation. However, instead of one ionization step 
being "frustrated" as is the case for pathway B, both ionization steps are "frustrated" in double 
H* formation. Moreover, we compute the screened nuclear charge that drives the explosion of the 
nuclei during double H* formation. 
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Many interesting physical phenomena arise during 
fragmentation of molecules driven by intense infrared 
laser fields. Examples of such phenomena are bond- 
softening and above-threshold ionization [l|, Q , molecu- 
lar non-sequential double-ionization and enhanced 
ionization [3|, l7l-[lQj. Another interesting phenomenon, 
the formation of single highly excited neutral fragments, 
has been recently observed in strongly driven H2 [ll| 
and other molecules p!2l-[l5|. The formation of excited 
fragments has been attributed to "frustrated" tunnel- 
ionization [l6|. In a very recent theoretical study of 
strongly driven H2 [l^]^ two main routes were identified 
leading to single H* formation. These two routes differ 
on whether the first (pathway B) or second (pathway A) 
ionization step is "frustrated" . 

In this work, using the theoretical model developed in 
[13], we study the formation of two highly excited H* 
atoms during the breakup of H2 by intense infrared laser 
fields. Multiple electron capture in Rydberg states was 
recently observed in strongly-driven Ar2 dimers [18]; no 
break-up channel with two neutrals as end fragments was 
reported in this latter study. We show that, as for path- 
way B of single H* formation, two-electron effects are im- 
portant in the pathway leading to double H* formation. 
Moreover, we quantify how the two highly excited elec- 
trons in double H* formation screen the nuclear charge 
and compute the effective charge that drives the Coulomb 
explosion of the nuclei. 

Our quasiclassical model was discussed in detail in ref- 
erence [13. It fully accounts for both the electronic and 
nuclear motion in order to accurately describe, among 
other processes, the formation of one and two H* atoms. 
Treating both electronic and nuclear motion in theoret- 
ical studies of strongly-driven molecules is a challenging 
task. This is the reason why, with few exceptions [19y,i,2Qj, 
in most studies either the nuclei are fixed [gI, [21| and only 
the electronic motion is considered or the electronic con- 
tinuum is ignored and only the nuclear motion is studied 

0. 

For completeness, we briefly describe our quasiclas- 



sical model. First, we set up the initial electronic 
phase space distribution. We consider parallel align- 
ment between the molecular axis and the laser elec- 
tric field (along the z axis) to complement our stud- 
ies of single H* formation in [17]. The field is taken 
to be E(t) = Eo(t) cos(cjt) at 800 nm corresponding to 
uo = 0.057 a.u. (a.u. - atomic units). In our simula- 
tion, the pulse envelope Eo(t) is defined as Eo(t) = Eq 
for < t < lOT and Eo(t) = Eocos2(cj(t - 10T)/8) for 
lOT < t < 12T with T the period of the field. We start 
the time propagation at cjto = (/>o where the phase of the 
laser field 0o is chosen randomly, see |23l-[25|. If the in- 
stantaneous field strength at phase (j)Q is smaller than 
the threshold field strength for over-the-barrier ioniza- 
tion, we assume one electron (electron 1) tunnel- ionizes, 
i.e., tunnels through the field- lowered Coulomb potential 
to the continuum whereas the initially bound electron 
(electron 2) is described by a one-electron microcanonical 
distribution. If the instantaneous field strength at phase 
00 allows for over-the-barrier ionization we use a dou- 
ble electron microcanonical distribution (see [25|). For 
both intensity regimes we use the tunneling rate pro- 
vided by the semiclassical formula in reference [26| with 
field strength equal to the instantaneous one at We 
use 0.57 a.u. (1.28 a.u.) as the first (second) ioniza- 
tion potentials. Second, for the initial phase space dis- 
tribution of the nuclei we use the Wigner function of the 
ground state (energy 0.01 a.u.) of the Morse potential 
[23; the parameters in the latter potential best describe 
the ground vibrational state of H2, for details see [13] • 

Third, we transform to a new system of "regularized" 
coordinates [2^, [29| . This transformation explicitly elim- 
inates the Coulomb singularity [25|. We then propagate 
the full four body Hamiltonian using the classical trajec- 
tory Monte Carlo method [30]. During the time prop- 
agation we allow the initially bound electron to tunnel 
at the classical turning points along the molecular axis 
using the WKB approximation [31, 32]. Finally we select 
only those trajectories that during the breakup of II2 re- 
sult in the H*+H* channel (where * denotes an electron 
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FIG. 1. Top row: The distribution of the quantum number n 
of the final Rydberg state occupied by either electron in the 
double H* formation channel. Bottom row: The final energy 
distribution of either of the H* fragments. The left column 
is for an intensity 1.5 x 10^^ W/cm^ and the right one is for 
2.5 X 10^^ W/cm^ 



300 
150 


-150 
-300 





H* 

e ^ — 5 


^e \ 

vny 


(b) Pathway B a/^ 




H2 




H2 a/^-'^'^' 


















\e J 




Ve J 





3 6 9 12 15 18 


21 


3 6 9 12 15 18 


21 




Time (field cycles) 




Time (field cycles) 





FIG. 2. (Colour online) Schematic illustration of the route 
leading to double H* formation (a) and that of pathway B of 
single H* formation (b). Shown is the time-dependent posi- 
tion along the laser field for electrons (black lines) and ions 
(grey broken lines). 
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in a n > 1 quantum state). To identify the electrons cap- 
tured in a Rydberg n quantum state of H* we follow the 
method described in reference (33| . 

In order to study the intensity dependence of dou- 
ble H* formation we have considered an intensity of 
1.5 X 10^^ W/cm^ in the tunneling regime and an inten- 
sity of 2.5 X 10^^ W/cm^ in the over-t he-barrier regime; 
however, for the latter intensity, most of the trajecto- 
ries are initiated with the tunneling model. We compute 
the final n quantum number [see Figure [TJa) and (b)] for 
either H* fragment for both intensities; at least 20000 
double H* events are considered. The final n quantum 
number distribution peaks at n = 8 for both intensities 
resembling the n distribution for singly excited neutrals 
formed either in strongly driven H2 [13 oi" strongly driven 
atoms [l6|. Moreover, the final energy distribution of ei- 
ther H* fragment in the H*-hH* channel [see Figure [TJc) 
and (d)] resembles that for the or H* fragments in 
the H*+H+ + e~ channel [l^- For the lower intensity 
[Figure ffl^c)] the energy distribution peaks at approxi- 
mately 3.5 eV. For the higher intensity [Figure [Tfd)] the 
peak shifts to higher energies since the nuclei Coulomb 
explode at smaller internuclear distances. 

In Figure [21(a) we show the pathway leading to dou- 
ble H* formation where electron 1 tunnel-ionizes very 
quickly, quivering in the laser field, while electron 2 
tunnel-ionizes after a few periods of the laser field. How- 
ever, when the field is turned off, both electrons 1 and 2 
do not have enough drift energy to escape and occupy 
instead a Rydberg state of an H atom. Similarly, in 
pathway B for single H* formation [l^, see Figure [2(b), 
electron 1 tunnel-ionizes very quickly, quivering in the 



FIG. 3. The distribution of the field phase 00 at the time when 
electron 1 tunnel-ionizes in the initial state for the double 
H* formation channel (•), and for pathway B of the single 
H* formation channel ( ). Figure (a) is for an intensity of 
1.5 X lO^^W/cm^ and (b) for 2.5 x lO^^W/cm^ 



field, while electron 2 tunnel-ionizes and escapes after a 
few periods of the laser field; however, when the field is 
turned off, electron 1 does not have enough drift energy 
to escape and occupies a Rydberg state of an H atom. 
Hence, while in the pathway leading to double H* forma- 
tion both ionization steps are "frustrated" in pathway B 
only the first one is. In what follows, we show that in all 
other respects the break-up channel leading to two H* 
atoms and pathway B of single H* formation are similar. 

We explore the route to double H* formation by, first, 
plotting in Figure [3] the distribution of the field phase, 
^0, when electron 1 tunnel- ionizes in the initial state. 
Two H* atoms form when ^0 is around 0° for the lower 
intensity [Figure [31(a)] and around -30° for the higher 
intensity [Figure [S^b)]. The reason double H* formation 
shifts from small values (extrema of the field) to larger 
values with increasing intensity is the onset of saturated 
ionization of the neutral molecule [25[. Comparing the 
distributions of the field phase (j)Q for pathway B of single 
H* formation and for double H* formation in Figure [31(a) 
(low intensity) and (b) (high intensity), we find that the 
two distributions are almost identical. 

Next, we ask how electron 2 gains energy to transition 
from the ground state of the H2 molecule to a high Ry- 
dberg state of an H-atom. We find that electron 2 gains 
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FIG. 4. (Color online) The 2-d distribution of the inter- 
electronic distance |fi — f2| versus the inter- nuclear distance 
R, at the time when the initially bound electron 2 tunnel- 
ionizes, for an intensity of 1.5 x 10^^ W/cm^ (a) and for 
2.5 X 10^^ W/cm^ (b). 
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FIG. 5. The mean inter-electronic distance (|fi — 12 
function of time, for an intensity of 1.5 x 10^^ W/cm 
double H*formation (solid), and for pathway B of single H 
formation (dashed line). 



energy through a strong interaction with the laser field 
that resembles "frustrated" enhanced ionization in 
("frustrated" since electron 2 occupies a Rydberg state 
instead of escaping). This resemblance is corroborated by 
i) the final energy distribution being similar for double 
H* formation [Figure [TJc) and (d)] and enhanced ion- 
ization [195 and ii) by our finding that electron 2 pref- 
erentially tunnel-ionizes when the nuclei are roughly 5 
a.u. apart as is the case for enhanced ionization. In 
Figure [H we plot the double differential probability of 
the inter-electronic distance versus the inter-nuclear dis- 
tance at the time when electron 2 tunnel-ionizes. Indeed, 
we find that electron 2 preferentially tunnel-ionizes when 
the nuclei are roughly 5 a.u. apart, with the inter- nuclear 
distance shifting to smaller values for the higher inten- 
sity. This is similar to our finding in [l^, for single H* 
formation. 

The questions that naturally arise next are whether 
electronic correlations are present while two H* atoms 
form and whether these correlations are similar in 
strength to the electronic correlation in pathway B of 
single H* formation. In the latter case, we have shown 
that electron-electron correlation mostly resembles the 
delayed pathway of NSDI Q (also referred to as re- 
collision-induced excitation with subsequent field ioniza- 
tion, RESI [111) where the electronic correlation is weak. 
For the delayed NSDI pathway, the re-colliding electron 
returns to the core close to a zero of the field, trans- 
fers energy to the second electron, and one electron es- 
capes with a delay more than a quarter laser cycle after 
re-collision. Comparing, in Figure O the mean inter- 
electronic distance as a function of time for the pathway 
where two H* atoms form and for pathway B of single 
H* formation, we find that the two mean distances are 
very similar. Indeed, a soft re-collision is present in both 
cases. Moreover, the two electrons stay closer to each 
other during double H* formation rather than single H* 
formation; this is to be expected since in the former case 
it is both electrons that undergo "frustrated" ionization 
oscillating in the vicinity of the nucleus. Thus, electron- 
electron correlation is similar in double H* formation and 



pathway B of single H* formation if not slightly stronger 
in the former case. 

Further, we find that the probability (out of all trajec- 
tories) for double H* formation reduces from 0.3% for the 
lower intensity to 0.2% for the higher one. Thus, forming 
two H* atoms is roughly 40-50 times more rare than form- 
ing one H* atom. The reduction of the probability from 
0.3% to 0.2% for double H* formation is consistent with 
a decrease, for increasing intensity, of electronic correla- 
tion in the form of re-collisions. Indeed, when comparing 
Figure [41(a) with (b) we find that at the time when elec- 
tron 2 tunnel-ionizes the two electrons are much closer 
to each other (smaller inter-electronic distance) for the 
lower intensity. 

We now gain further insight into the dynamics of the 
Coulomb exploding nuclei by plotting in Figure [6] the 
mean total kinetic energy of the nuclei (Ekin) versus the 
mean inverse inter- nuclear distance (1/R). Note that in 
Figure[6]time increases as (1/R) decreases. From early on 
up until the time (roughly 3.8 laser cycles) when electron 
2 tunnel-ionizes at (1/R) ^ 0.15 a.u. it is mostly electron 
2 that significantly screens the Coulomb repulsion of the 
two nuclei. This screening is corroborated by the small 
slope of the (Ekin) versus (1/R) curve for distances up to 
(1/R) ~ 0.15 a.u. After electron 2 tunnel-ionizes the two 
nuclei move fast away from each other. Indeed, the slope 
of the (Ekin) versus (1/R) curve increases for (1/R) > 
0.15 a.u to 0.86 for both intensities. 

This slope of 0.86 quantifies how much the two elec- 
trons that remain bound to highly excited states screen 
the Coulomb repulsion of the exploding nuclei at large 
inter- nuclear distances. To show that this is the case, 
we plot (Ekin) versus (1/R) for break-up channels where 
both electrons escape. Such channels are the double 
ionization through re-collision or through enhanced ion- 
ization with two H+ ions and two escaping electrons as 
the final fragments. To identify double ionization events 
through re-scattering, we use as a rough criterion the con- 
dition that electron 2 escapes without tunnel-ionizing. If, 
however, electron 2 tunnel- ionizes, we register these tra- 
jectories as double ionization through enhanced ioniza- 
tion events. Figure [6] shows that the two nuclei move for 
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scattering rather than through the enhanced ionization 
channel. This is also expected since in the former case 
electron 2 escapes at times smaller than the time it takes 
for electron 2 to tunnel-ionize in the latter channel. For 
large times (small (1/R)), after the laser pulse is off, we 
find that the slope of the (Ekin) versus (1/R) curve for 
the double ionization channels is 1, consistent with two 
H+ nuclei Coulomb exploding unhindered by electronic 
motion. 



FIG. 6. The mean total kinetic energy of the nuclei (Ekin) 
versus the mean inverse inter- nuclear distance (1/R), for 
1.5 X 10^^ W/cm^ (a) and for 2.5 x 10^^ W/cm^ (b). The 
solid line is for the double H* formation channel, the dashed 
black/gray line is for the fragments in the double ioniza- 
tion channel through re-scattering/enhanced ionization. The 
arrows indicate the (1/R) corresponding to a time equal to 
12T (duration of the laser pulse). The insets show the seg- 
ments we fit for times larger than 12T to obtain slopes of 0.86 
for the double H* formation channel and of 1 for the double 
ionization channels. 



all times faster away from each other in the double ion- 
ization channels rather than in the double H* formation 
channel. This is expected since in the former case both 
electrons escape fast, not screening the Coulomb repul- 
sion of the nuclei. Moreover, the nuclei move faster away 
from each other in the double electron escape through re- 



Concluding, we have found that in double H* forma- 
tion both tunnel- ionization steps are "frustrated". We 
have shown that the route to forming two H* atoms is 
a pathway similar to pathway B of single H* formation. 
Specifically, two-electron effects are present in both sin- 
gle and double H* formation channels. That is, in dou- 
ble H* formation, electron 2 gains energy through a weak 
interaction with electron 1 resembling delayed NSDI in 
H2. Moreover, electron 2 gains mostly energy through a 
strong interaction with the laser field resembling "frus- 
trated" enhanced ionization of H2", while electron 1 occu- 
pies a high Rydberg state. Two-electron effects diminish 
with increasing intensity with electron 2 gaining energy 
mainly through an interaction with the laser field. We 
emphasize that our 3-d quasiclassical formulation for de- 
scribing break-up channels during Coulomb explosion of 
strongly driven H2 is general and can be applied to study 
the break-up of driven multi-center molecules. 
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